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Using a recently proposed gauge invariant formulation of light-cone averaging, together with 
adapted "geodesic light-cone" coordinates, we show how an "induced backreaction" effect emerges, 
in general, from correlated fluctuations in the luminosity distance and covariant integration measure. 
Considering a realistic stochastic spectrum of inhomogeneities of primordial (inflationary) origin we 
find that both the induced backreaction on the luminosity-redshift relation and the dispersion are 
larger than naively expected. On the other hand the former, at least to leading order and in the 
linear perturbative regime, cannot account by itself for the observed effects of dark energy at large- 
redshifts. A full second-order calculation, or even better a reliable estimate of contributions from 
the non-linear regime, appears to be necessary before firm conclusions on the correct interpretation 
of the data can be drawn. 

PACS numbers: 98.80-k, 95.36.-fx, 98.80.Es 

I. INTRODUCTION 

The so-called concordance (or ACDM) model, based on a suitable combination of dark matter, dark energy and 
baryons for an overall critical density, has become the reference paradigm for the late - i.e. post-equality epoch - 
evolution of our Universe (see e.g. [ ]). It accounts equally well for the CMB data, the Large Scale Structure and, 
even more significantly, for the supernovae data in terms of a cosmic acceleration [2] . 

Strictly speaking these three tests of the concordance model are not at the same level of theoretical rigor. While 
the first two have to do, by definition, with the inhomogeneities present in our Universe, the third is based on an 
ideal homogeneous and isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) geometry. It is clear that a better 
treatment of cosmic acceleration should take inhomogeneities into account, at least in an average statistical sense. Only 
when this is done we can establish in a convincing way whether ACDM gives a simultaneous consistent description of 
the above-mentioned body of cosmological data. 

This realization has led to a vast literature about averaging cosmological observables in realistic inhomogeneous 
cosmologies (see e.g. [■>] for recent reviews). The conclusions, however, are still rather controversial: according to 
some authors [4] present inhomogeneities might explain, by themselves, cosmic acceleration without any need for 
dark-energy contributions; according to others [•">] the effect of inhomogeneities is, instead, completely negligible. The 
truth may lie somewhere in between, in the sense that a quantitative understanding of inhomogeneities effects could 
be important in order to put precise constraints on dark-energy parameters, such as the critical fraction of dark-energy 
density, Qa, and the time evolution of its effective equation of state, wa{z). 

In the first papers studying the dynamical effects of averaging, the problem was approached mainly following 
Buchert's prescriptions [>'], namely averaging inhomogeneities over spacelike hypersurfaces and computing the ensuing 
"backreaction" on the averaged geometry. Nonetheless it is clear that a proper treatment of cosmic acceleration - 
which is indeed revealed through the experimental study of the luminosity-distance to redshift relation [2] - needs 
to consider the backreaction of averaged inhomogeneities along our past light cone, which is a null three-dimensional 
hypersurface. 

While a number of papers have emphasized the importance of light-cone averaging (see e.g. [7-!i]), a fully covariant 
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and gauge invariant formulation of such a procedure was only given recently [10], by generalizing to null hypersurfaces 
an analogous prescription previously derived for spacelike hypersurfaces [II, 12]. Actually, it turns out that the 
physically meaningful (covariant and gauge invariant) average of a scalar over a null hypersurface reduces to averaging 
over an appropriate two-dimensional surface a scalar object which is non local, as the integrand itself contains integrals 
along lightlike geodesic curves lying on the given null hypersurface. Nonetheless, we shall simply refer to this procedure 
as "hght-cone averaging" [10]. 

The aim of this paper is to apply such an averaging procedure to the luminosity distance of a light source lying 
on our past light-cone, in order to compute the possible backreaction of inhomogeneities on the luminosity-redshift 
relation. We will show the emergence of an effect - called hereafter "induced backreaction" - which arises from a generic 
correlation between the inhomogeneities present in the variable we want to average (e.g. the luminosity distance) and 
those appearing in the covariant integration measure"'^ . We stress immediately that our induced backreaction accounts 
only for a part of the total effect of inhomogeneities. As we shall discuss in detail below, the leading-order induced 
backreaction can be computed in terms of linear perturbation theory, while a complete leading-order calculation (left 
to future work) would require perturbation theory up to second order. 

We will use a simple phenomenological model of inhomogeneous geometry based on a spatially flat FLRW metric 
which includes, to first order, scalar perturbations of primordial (inflationary) origin. Such perturbations are often 
conveniently parametrized in the longitudinal (or Newtonian) gauge [13]. However, we will take advantage of the 
gauge invariance of our formalism to use an adapted system of coordinates defining the so-called "geodesic light-cone" 
(GLC) gauge introduced in [loj. In that gauge the light-cone averages of luminosity distance and redshift take very 
simple, exact expressions that keep all the required degrees of freedom for being applicable to general geometries. 

By further taking the ensemble average of the inhomogeneous terms, and using the stochastic properties of our 
model of perturbations, we will compute to leading order the induced backreaction on the luminosity distance, and 
the corresponding variance. If we limit ourselves to first order computations in the range of scales where linear 
perturbation theory is reliable, we find that the induced (second-order) corrections to the luminosity-redshift relation 
of the homogeneous CDM model, although much bigger than one could have naively expected, are not large enough 
to mimic a sizable fraction of J^a- However, our formalism suggests that other second-order corrections could be 
even larger, thus confirming the importance of performing a full second-order calculation or, even better, a reliable 
calculation in the non-linear (short scale) regime. 

The paper is organized as follows. In Sect. II we recall definition and basic properties of GLC coordinates and 
give exact, non-perturbative expressions for the redshift, the luminosity distance and the light-cone average of the 
latter at constant redshift. In Sect. HI we give the linear order transformations expressing the metric of a perturbed 
FLRW geometry in GLC form, starting from the more commonly used Newtonian gauge. We then express the 
luminosity distance as a function of the redshift and of the angular position of the source, to first order in the given 
perturbed geometry. In Sect. IV we take the light-cone average of the computed luminosity distance and we show 
how backreaction effects automatically emerge, in general, from correlations between the averaged variable and the 
covariant integration measure. We explicitly define such an induced backreaction and the cosmic variance, to leading 
order, taking into account the stochastic properties of our model of perturbations and thus including the appropriate 
ensemble average of the inhomogeneous corrections. In Sect. V we implement analytical calculations of the induced 
backreaction terms and of the dispersion, and in Sect. VI we present the corresponding numerical results for a specific 
spectrum of primordial scalar perturbations and for a given model of transfer function. Our main conclusions are 
finally summarized in Sect. VII. We also present, in Appendix A, a sample of the analytic manipulations needed for 
reducing our (light-cone plus stochastic) averages to the form of a one-dimensional integral over the wavenumber k 
labelling the scalar perturbation spectrum. 

II. THE GLC GAUGE AND THE AVERAGED LUMINOSITY-REDSHIFT RELATION 

A. Reminder of the GLC gauge and of its main properties 

Our purpose, in this paper, is to compute the light-cone average of the luminosity distance as a function of the 
redshift by applying the gauge invariant procedure introduced in [10]. One of the main virtues of using a gauge 
invariant formalism is the freedom of choosing a gauge particularly adapted to the problem at hand. In the case 
of spacelike averaging, for instance, a convenient coordinate system corresponds to a gauge where the averaging 



As discussed in [Id], this integration measure will also induce "backreaction" terms of the type usually discussed in the literature [R], 
namely terms that arise from (generalized) commutation rules between differential operators and averaging integrals. 
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hypersurfaces are identified with constant-time hypersurfaces. In many appUcations the averaging hypersurfaces are 
chosen indeed as the ones associated to a class of geodesic observers corresponding to constant values of the time 
parameter t of the synchronous gauge (see e.g. [14]). 

Similarly, for light-cone averages, it is convenient to identify the null hypersurfaces with those on which a null 
coordinate takes constant values. For this reason we have introduced in [Lii] an adapted system of coordinates 
- defining what we have called a "geodesic light-cone" (GLC) gauge - where the averaging prescription greatly 
simplifies, while keeping all the required degrees of freedom for applications to general geometries. Furthermore we 
are also able to identify the timelike coordinate of the GLC gauge with the cosmic time t of the synchronous gauge. 
As a consequence we can easily introduce a family of geodesic reference observers which exactly coincide with the 
static ones of the synchronous gauge. 

Let us first recall, as discussed in [ ], that the coordinates a;^ — {t,w,9"-), a — 1,2, specifying the metric in the 
GLC gauge correspond to a complete gauge fixing of the so-called observational coordinates, defined e.g in [7, 15, 16]. 
The GLC metric depends on six arbitrary functions (T, a two-dimensional "vector" and a symmetric matrix jab), 
and its line element takes the form 

ds^i^c = '^^dw^ - 2TdwdT + jabidr - U''dw){d9'' - U^dw). (2.1) 
In matrix form, the metric and its inverse are then given by: 

/ -1 -T-i -?7VT\ 
- -T-^ , (2.2) 

V-(C/")^/T 0^ T"^" ) 

where = (0,0), Ub = (J7i,[/2), while the 2x2 matrices 7ah and 7"* lower and raise the two-dimensional indices. 
Clearly w is a null coordinate (i.e. d^wd^w ~ 0), and a past light-cone hypersurface is specified by the condition w — 
const. We can also easily check that S^r defines a geodesic flow, i.e. that (S^t) V ^ (9^t) = (as a consequence of 
the relation g'^'^ = — 1). 

In the limiting case of a spatially flat FLRW geometry, with scale factor a, cosmic time i, and conformal time 
parameter 77 such that drj = dt/a, the transformations to the GLC coordinates and the meaning of the new metric 
components are easily found as follows ["']: 

T — t, w — r + ri, T = a(t), 

U" = 0, labdO^-de^ = a^{ty{d0^ + sin^ Odcj)^). (2.3) 

Even though we will be mainly using the GLC gauge for a perturbed FLRW metric in the Newtonian gauge, it is 
important to stress that the equality between the coordinate t and the proper time t of the synchronous gauge holds 
at the exact, non perturbative level: it is always possible, in fact, to choose the GLC coordinates in such a way that 
T and t are identified like in the above FLRW limit. 

In order to illustrate this point let us consider an arbitrary space-time metric written in the synchronous gauge 
(with coordinates — (i,X*), i,j ~ 1,2,3), where the line element takes the form: 

dslc = -dt^ + h.jdX'dX^. (2.4) 

Let us impose the condition t = t, and check whether we run into any contradiction with the exact metric transfor- 
mation 

Using Eq. (2.2) for the GLC metric, and considering the transformation for the g^^Q components, we then obtain the 
conditions 

gfa = {-1.0} = -[dr + r-\d.^ + W'da)] - u'^d.x^^ = (2.6) 

where = —d^r is the four- velocity of the geodesic GLC observer, and where A denotes an affine parameter along the 
observer world-line. So the requirement t = t boils down to the statement that along the geodesic flow of the vector 
field the SG coordinates X' are constant. This clearly defines the coordinate transformation in a non-perturbative 
way, and also shows that the geodesic observer = —d^r of the GLC gauge corresponds to a static (and geodesic) 
observer in the synchronous gauge. It follows that the identification t = t can always be taken for any space-time 
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metric, and that this simple connection between GLC and synchronous gauge has vahdity far beyond the particular 
FLRW case or its perturbed generalizations. 

We also remark that, in GLC coordinates, the null geodesies connecting sources and observer are characterized by 
the simple tangent vector /c^ = g^'^d^w = g^"^ = — J^T"^, which means that photons reach the observer travelling at 
constant w and 0°" . This makes the calculation of the area distance and of the redshift particularly easy. Consider, 
for instance, a light ray emitted by a static geodesic source at the intersection between the past light-cone of our 
observer, w = wq, and the spatial hypersurface t = r^, and received by such static geodesic observer at r = tq > t^. 
The associated redshift Zg is then given by [ ]: 

a + z) = (^""^^^ ^ idf^^d^r), ^ TK,To,r) 

where the subscripts "o" and "s" denote, respectively, a quantity evaluated at the observer and source space-time 
position. The expression for the angular distance will be explicitly derived in Sect. II B. 

Let us finally recall that, in GLC coordinates, the covariant average of a scalar quantity S{t, w, 9°") over the compact 
two-dimensional surface S, defined by the intersection of our past light-cone w — Wq with the spacelike hypersurface 
T = Tg, is simply given by [10]: 

^ d^x^5{w - wa)6{T - Ts)S{t, w, r ) \d^Td^'w\ 

= / ' (2-8) 

where 7 — det7a6. In the case of interest for this paper, namely light-cone averages on surfaces of constant redshift 
z — Zs, one then obtains [10] 

_ / d^0^j{wo,T{z,,wo,9-),e^~)S{wo,T{zs,w„,0-),9'') 
/ d^9^j{wo,T{zs,wo,9-),9b) 

where t[zs,wq,9°') has to be determined by solving the redshift equation (2.7) for as a function of wq, tq, Zg and 
9°'. This general result will now be applied to the case in which S is identified with the luminosity distance dL- 



B. Light-cone average of the luminosity distance 



Let us first recall that the luminosity distance c?^ of a source at redshift z is related in general to the angular 
distance dA of the source (as seen from the observer) by the so-called Etherington (or reciprocity) law [17]: 



dL = [l + zfdA ■ 



(2.10) 



In the particular case of an unperturbed, spatially flat FLRW background, and for a source with redshift Zs, the 
distance dA is simply given by 



jFLBW 

"A 



(2.11) 



where = a{r]s), while tjq — rj^ denotes the conformal time interval between the emission and observation of the light 
signal. For the unperturbed metric, on the other hand, we have 1 + z = ao/a{t), and drj = dt/a = —aQ^dz/H, where 
H = d{\na)/dt. Hence: 



^FLRW 



Vo 



(zs) = (1 + Zs)ao / d?7 = (1 + Zs) 



dz 
H{z) 



Ho 



dz 



3(H 



(2.12) 



In the last equality we have used the standard (spatially flat) Fricdmann equation for H (see e.g. [ ]), assuming 
that the given homogeneous model has perfect fluid sources with present fractions of the critical density Qno Sind 
barotropic parameters Wn- Expanding in the limit — >■ we also obtain the expression 



^FLRW 



(Zs) 



1 



Ml- 3^ w„ri„, 



z^ 



0{zl 



1 

Ho 



Zs + \{l-qo)zl + 0{zl) 



(2.13) 
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which shows the well known sensitivity of the term quadratic in Zs to the composition of the cosmic fluid through the 
deceleration parameter go- 

Let us now discuss how this well known result for is modified when including generic inhomogeneities. We recall, 
to this purpose, that in a generic metric background the angular distance c?a can be computed by considering the null 
vector A;^ — dx^ /d\ tangent to the null ray connecting source and observer (i.e. belonging to the congruence of null 
geodesies forming the observer's past light-cone). Here A is an affine parameter along the ray trajectory, chosen in 
such a way that A = at the observer and A = A^ at the source position. The expansion Q of the congruence of null 
rays is then given by G = V^fc^, and the corresponding angular distance dA is defined by the differential equation 
[15, 19]: 



^il.dA)=^ = ' 
d\^ ' 2 2 



{IndA) = - = -V^¥\ (2.14) 



Consider now a generic metric in the GLC gauge, and the null vector fc^ = d^w such that fc^ — dx^'/dA = — (5^T. 
We have, in this case, 

V^k'^ = ^T~'dr (In V7) = ^ (In Vl) , (2-15) 
where 7 — det7ab, and the integration of Eq. (2.14) gives, in general. 



4(A) = c/W), (2.16) 

where c is independent of A. In order to fix the constant c we may recall the boundary conditions required to guarantee 
that the null rays generated by belong to the past light cones centered on the world line of our observer. In the 
limit — >■ (or A 0), where 9'^ 9°- (see in particular next section, Eq. (3.8)), such conditions require in 
particular that [I''] 

lim ^ = sine}, = sin^\ (2.17) 

A-).0 d\ 

where, in the last equality, we have used the already mentioned fact that 6*" is constant along the null geodesic. This 
clearly fixes c = 1/ sin 6*^, and uniquely determines the angular distance as 

d^(A) =7i/4(A) (sin^i)"'^'. (2.18) 

Inserting this result into Eq. (2.9) we finally arrive at an exact expression for the light-cone average of the luminosity 
distance, as a function of z^, in the GLC gauge: 

, ,2 1 d'O I'/HwQ, Tjzs, wo, g°), e')dA{wo, r(z„ ^0, e-),9'') 

{aL)wo,z, - (i- + Zs) „ ,„~ , ; ~ , ~, , ■ l^-J-yj 

/ d^e j^/^wo,T{z„wo,9^),eb) 

This equation, being exact, can be applied in principle to any given highly inhomogeneous cosmology. Specifically, 
it can be applied to a Lemaitre-Tolman-Bondi (LTB) model [-. ] even in the case of a generic observer shifted away 
from the symmetry centre of the isotropic geometry. In this paper, however, we shall be mainly interested in working 
out the consequences of light-cone averaging for a particular model of perturbed FLRW Universe. Since we know 
how to describe the latter, for example, in the Newtonian (longitudinal) gauge, we need in general to connect such a 
coordinate system to the GLC gauge. 

III. LUMINOSITY DISTANCE IN A PERTURBED FLRW GEOMETRY 
A. First-order coordinate transformation from the Newtonian to the GLC gauge 

Let us consider the scalar perturbations of a conformally fiat FLRW background to describe the particular model 
of inhomogeneous geometry we are interested in this paper. Choosing in particular the so-called longitudinal or 
(conformally) Newtonian gauge (NG), using spherical coordinates (r, 9°") = (r, 9, (f>), and going up to the first order in 
perturbation theory, it is well known that the model is parametrized by the following (inverse) metric tensor [I'!]: 

9%G = a^'iv) diag (-1 + 2$, 1 + 2*, (1 + 2^^)70"^) . (3.1) 
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Here 



^ diag {r~^ ,r'^ sin-^ 9) 



(3.2) 



and $, ^' are the usual gauge invariant Bardeen potentials [ I >>], general functions of 77, r, 9 and 0. We will assume, for 
simplicity, that the matter sources have vanishing (or negligible) anisotropic stress, so that (f> = 5*. 

For the subsequent computations - in particular, for the application of the simple average prescription given 
previously - we need to re-express this metric in GLC form, finding the transformations from the NG coordinates 
yf^ — {ri,r,9,(p) to the GLC coordinates x'^ — {t,w,9^,9'^), and computing the reparametrized metric as 



drP Br" 

To this purpose we will introduce the useful (zeroth-order) light-cone variables r]± = 77 ± r, such that 



(3.3) 



V ^ 2 {11+ + V- ) 
with corresponding partial derivatives 



(3.4) 



(3.5) 



Using these variables we solve the three differential equations obtained from Eq. (3.3) for the components goj^c = 
— 1, Qq^q = 0, gQ^Q = 0, by imposing the boundary conditions that i) the transformation is non singular around the 
observer position at r = 0, and ii) that the two-dimensional spatial sections r = const are locally parametrized at 
the observer's position by standard spherical coordinates, i.e. 9"'{0) = 9"" = {9^(f)). The sought for transformation can 
then be written, to first order in ^P, as follows. 



11+ + / (ia;*(77+,x,r) , 

ah I 



(3.6) 
(3.7) 
(3.8) 



v+ 



where ^{ri+,r]-,9°') = ^{r],r,9'^),^'^^{ri+,ri-,9°') = ^"'^{rj,r,9°') and ?7i„ represents an early enough time when the 
perturbation (or better the integrand) was negligible. We can easily check that, to zeroth order in ^E*, we recover the 
homogeneous transformation (2.3) as expected. 

To first order in ^E" we can use again Eq. (3.3) to compute the non-trivial entries of the GLC metric (2.2), and 
obtain: 



1 



l + *(77+,77+,r)- / dxd+^{ii+,x,9'') 

1 



+ / d77'a(7/)9.*(77',r,r); 



(3.9) 



= / dxdb^{v+,x,9'') - —70"" / d7j'a{rj')dt'f{v',r,9'') 



dx 



v+ 



a{r]) 



jS%V+,x,0^) I dy dt^irj+,y,9'^) 



- lim 

2 2;— i.r;+ 



jr{v+,x,9^) I dy dt^iv+,y,9'^) 



[l + 2vI/(,7,r,r)]7o"' 



1 



dxd, ^l\ii+,x,9"') / dy9d^(r;+,2/,r) +ao6 



■ (3.10) 
.(3.11) 



The term '5(77+, 77+, 0°), appearing in the general equation for T, denotes (at any rf) the value of the perturbation 
potential evaluated at the tip of the light-cone connecting the origin (r = 0) to the point y^^ = {ri,r,9'^) (namely, 
^P(r,+ ,7?+,r) = vI,(^ + r-,0,r)). 
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The above transformations can be immediately applied to obtain an explicit expression for the redshift parameter 
Zs- By inserting the result (3.9) into Eq. (2.7) and considering that the source is located on the past light-cone of the 
observer (identified by the equation w = wq = ?7o)i we obtain: 



1 



1 + j(2„r) 



where J — Ix 



Ir, and where: 

/+ = 



I. 



"^^ a(n') 



a('7s) 



a(n') 



(3.12) 

(3.13) 
(3.14) 



We have defined rj^ = rjs ir^, — '^{Vs, Vo ~ Vs, ^°)i '^o = '^{vo, 0, ^'^) and we have used the zeroth-order light-cone 
condition i]^ = V+ = Vo- It should be stressed, however, that while the integrals appearing in Ir are evaluated at 
constant r (namely along timelike geodesies), all the other integrals are evaluated at fixed rj^ (i.e. along null geodesies 
on the observer's past light-cone). 

The contribution associated to Ir can also be rewritten as 



Ir = {Vs - Vo) ■ n, 

where n is the unit tangent vector along the null geodesic connecting source and observer, and where 



drj 



V*(77',r, 9") 



(3.15) 



(3.16) 



are the "peculiar velocities" of source and observer associated to a geodesic configuration perturbed up to first order 
in the NG gauge. In a realistic situation one should add to Eq.(3.12) similar terms taking into account a possible 
intrinsic (non-perturbative) motion of source and observer, unless our theoretical predictions are compared with data 
already corrected for these latter Doppler contributions. Note also that our first-order expression for the redshift 
(3.12), valid in general for any given scale factor a{r]), is in full agreement with the expression obtained in Eq.(38) of 
[21], where Zs is computed for the particular case of a CDM-dominated model. 



B. The luminosity distance including first-order scalar perturbations 



We now apply the above coordinate transformation to find, in the Newtonian gauge and to first order in perturbation 
theory, the other relevant quantities for this paper. Let us start with the determinant 7 appearing in the angular 
distance (2.18). For a source emitting light at time rjg and radial distance we obtain from Eq. (3.11), to first order^. 



7-i(A,) = det7'^''(A,) = (a^r^sin^)-^ 1 + 4^, + 4 J2(z„ ^ ) 



(3.17) 



where: 



~ 1 Z"*- 



(the latter equality follows upon a simple integration by parts). Hence, from (2.18): 

/ . ^1 \ -1/2 



J2 



sin^^ 



(3.18) 



(3.19) 



Note that, to first order, we have 7 ^ = (711722) ^, and that, for these diagonal matrix elements, the operator di, commutes with 7q' 
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The last factor can be easily computed, to first order, by using Eq. (3.8) and the fact that the 9°" are constant along 
the null geodesic. It is easy to check that it amounts to a redefinition of J2, and that the above angular distance 
becomes 



where: 
J2 



1 



drj 



dA{Xs) = asr s[l - '^s - J2izs,0'')] , 
9^ + cot ede + (sin 0) -2^21 ^(^'^ ,^0 - r;', ^ ) 



1 



'70 - Vs J,,, Vo-V 



(3.20) 



rf^IL^As^- (3.21) 



(here A2 is the two-dimensional Laplacian operator on the unit two-sphere). 

For the full explicit expression of the luminosity distance cIl at constant rcdshift what we need, at this point, is 
the first-order expansion of the factor QsTs = a{ris)rs appearing in Eq. (3.20). To this purpose we start from Eq. 
(3.12), considering Zg as a constant parameter localising the given light source on the past light-cone (w — Wq) of our 
observer, and we look for approximate solutions for rjs = rj^^ZsjO""). 



Let us first define the zero-order solution 77^*^^ through the exact relation 



a{r). 



(0)n 



1 



where ao = 0(770). Expanding (3.12) with respect to the parameter 67] = rjg — rji^'^we then find: 



1 



0(77, 



(0)^ 



1 + Zs ao 



-[1+Hs ^77- J(z„r)] 



1 



l + z, 



■[I + H.577- J(z,,r)] , 



(3.22) 



(3.23) 



where Tis ~ d{\na{ris ))/dr]s , so that: 

Hs577- J(z„r) . (3.24) 

On the other hand, by applying Eq. (3.7) to the light-cone w = wq at the source position, we readily obtain 

u;o=77;-2A77*av = %, (3.25) 

where we have introduced the zero-order quantity Ar/ = 770 — r]i^\ and denoted by ^'av the average value of ^E* along 
the (unperturbed) null geodesic connecting source and observer: 

dx ^(77;, X, 0°) = -2 / d77' ^-(7/, 7/0 - 7/, e'') = -2A7/*av. (3.26) 

Combining this result with Eq. (3.24) we can then determine, to first order, the value of the radial coordinate rg{zg, 9°-) 
corresponding to the given rcdshift z^: 



r,(z„ 0") = ?i;o - ?7i°^ (^.) - (577 + 2A77«'av = A77 



1 L^.! L -x- 0\\i 

T / A I ^ ^ av 



Proceeding in the same way for we obtain, from Eqs. (3.12), (3.22), 

a,(z„(?'^) =a(7,f )[1 + J(z„r)], 
which, in turn, allows us to compute the value of asVa on the constant-z^ 2-surface: 

1 



[a,r,](z„r)=a(77(0))A77 



1 + 2*, 



1 



H,A77 



J{zs.9"-) 



The angular distance (3.19), for a source at rcdshift z^, can now be written as'^ 

1 



dA(z„(?'^) - 0(77(0) )A77 



1 + 2*, 



1 



H,A77 



J~*(7;„77o-77„r)- J2 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



Let us note that in the first order terms we can always safely identify ris with its unperturbed value rj. 



(0) 
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The (first-order, non-homogeneous, non-averaged) expression of in our perturbed background, referred to the 
unperturbed value (2.12), is finally given by: 

- l-*„..,,„-,„n + 2*„,+ (l-^),/-.,. (3.31) 

If we apply this general result to the particular case of a CDM-dominated Universe we find almost full agreement 
with the result for the luminosity distance at constant redshift computed in [21]. After several manipulations, in fact, 
it turns out that our is equivalent to the one found in [21], modulo a term which can be written as vq • n (in the 
notations of Eq. (3.15)). Such a term gives a subleading contribution to the backreaction and can be neglected with 
no impact on our final results. 

It should be noted, however, that by expanding the above expression in the small Zg limit, and comparing the 
result with the analogous expansion of the homogeneous distance (2.13), we could easily introduce a redefined value 
of Hq, say Hq'^'^, such that its inverse corresponds to the coefficient of the linear term in Zg for the perturbed relation. 
With such a "renormalized" Hubble parameter we have that H^'^^dL tends smoothly to H^d^^^^ for Zg — )■ 0, and 
we recover full agreement with . ] provided the same renormalization is applied there. Such a renormalization of 
is also suggested by (and closely related to) the first-order computation of the scalar expansion factor V^u^ for the 
flow wordlines = d^T of a local geodesic observer. 

At this point we could go on by taking the average of Hg°'^dL{zs, 0°), computing the associated backreaction, and 
evaluating the corrections to the standard homogeneous relation given by Hodf^'"^^ (zs). We have performed that 
exercise, but we have found that the contribution of the renormalization terms give large Zg-independent contributions 
to the variance. Furthermore, we think that renormalizing Hq at = is physically incorrect since, at very small- Zg, 
the backreaction is dominated by short-scale inhomogeneities which are deeply inside the non-linear regime (where 
even the concept of a Hubble flow becomes inappropriate) . We could instead try to renormalize Hq at some small but 
finite Zg, e.g. at a redshift corresponding to the closest used supernovae (say Zg ~ 0.015, see e.g. [22]), but then the 
results (although much better behaved) would depend on the choice of the particular "renormalization point" . Thus, 
it seems best to consider just the full "unrenormalized" expression (3.31) in a limited region of Zg where one can trust 
the approximations made, and use that expression for a phenomenological parametrization of the backreaction effects 
which could possibly include a redefinition of Hq. 



IV. COMBINING SPACE-TIME AND ENSEMBLE AVERAGES 



In the following sections the inhomogeneous deviations from the standard FLRW quantities are sourced by a 
stochastic background of primordial perturbations, satisfying ^ = 0, ^'^ 0, where the bar denotes statistical (or 
ensemble) average (see Sec. V). Hence, if we limit ourselves to a first-order computation of d^, we immediately obtain 
= df^^^^ . Non-trivial effects can only be obtained from quadratic and higher-order perturbative corrections, or 
from the spectrum of the two-point correlation function dL{z,d'^)dL{z' ,6''^) (discussed in detail in [2)]). 

In this paper we will consider the ensemble average not of d^ but of (d^), where the angular brackets refer to the 
light-cone average defined in Eq. (2.19) (see e.g. [23-25] for previous attempts of combining ensemble average with 
averages over spacelike hypersurfaces) . We will see that the light-cone average automatically induces quadratic (and 
higher-order) backreaction terms, due to the inhomogeneities present both in d^ and in the covariant integration 
measure, and since the ensemble average of such terms is non- vanishing we obtain, in general, (d^) 7^ df^^^^ . 

We will start this section with some general considerations on how to combine space-time and ensemble averaging, 
and how to isolate those terms in (d^) that we may genuinely call "backreaction" effects, i.e. effects on averaged 
quantities due to inhomogeneities. We shall also discuss how to estimate the variance around mean values due to such 
fluctuations^. Many of these considerations can be certainly found elsewhere, but are nonetheless presented here for 
the sake of being self-contained. 

Let us consider a typical average over the compact surface S (topologically equivalent to a two-sphere) embedded 
on the past light-cone w = wq at constant Zg. We simply denote such an average by: 

iSh = ^ , (4.1) 



* The importance of the cosmic variance for a precise measurement of the cosmological parameters, taking into account backreaction 
effects from averaging on domains embedded in a spatial hypersurface (according to [o]), has been recently pointed out also in [-"]■ 
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where cPfi is the appropriate measure provided by our gauge invariant prescription (see Eq. (2.9)) and S is the 
(possibly non local) scalar observable (di in our case). We can conveniently extract, from both (P^ and from S, a 
zeroth-order homogeneous contribution by defining: 

dV= (rfV)^°ni + Ai), S^S^°\l + <j), (4.2) 

and use the possibility of rescaling both integrals in Eq. (4.1) by the same constant, in order to normalize / ((i^/i)(°) = 1. 
We then easily get: 

/ S \ ^ /(d^M)W(l + M)(l + a) ^ ((l + M)(l + ^))o 

\sio)/ /(d2^)(o)(l + ^) l + (^)„ ' ^■■'> 

where we have dropped, for simplicity, the subscript of the averaging region S, and where we have defined by (. . . )o 
averages with respect to the unperturbed measure {(PfiY^^^ (we shall drop the subscript hereafter). 

Let us now perform the ensemble average, denoted by an overbar, paying attention to the fact that ensemble 
averages do not factorize, i.e. AB A B. A simple calculation leads to: 



(5/5(0)) ^ 1 + ((^) + + . (4.4) 

This last equation is supposedly exact but, as such, pretty useless. It becomes an interesting equation, though, if we 
can expand the quantities /i and cr in a perturbative series: 

M = cr^^cTt, (4-5) 

i i 

and we further assume that the first order quantities /^i,(Ti have vanishing ensemble averages (as it is the case for 
typical cosmological perturbations coming from inflation). In that case we can easily expand the result and obtain, 
for instance: 



(S'/S'(o)) ^ 1 + (ct2) + IBR2 + (0-3) + IBR3 + . . . (4.6) 

where 

IBR2 = (Aii<7i> -(m)(ai>, (4.7) 



IBR3 = (/^20-l) - (a*2)(cti) + (Micr2) - (M1>('^2> - (/il)(/iicri) + (MiXMiXc^i), (4-8) 

and where we have used again the non-factorization property, i.e. {fii){ai) ^ (^i) (cti), and so on. We see that 
the result contains, to a given order, both terms that depend on expanding S to that order (but not on the precise 
averaging prescription), and "induced backreaction" (IBR) terms that depend on correlations between the fluctuations 
of S and those in the measure. These latter terms only depend on lower-order perturbations of S and the measure 
separately. In particular, our first-order calculation provides the full second-order IBR effect that comes from the 
above interplay of and cti, although the full second-order result needs also the harder computation of a2 (but not 
of fi2)- Note also that, whenever the fluctuations of S and (Pn are uncorrelated, all IBR effects drop out. We will see 
below how to apply the above general reasoning to the particular case of di. 

Let us now discuss instead the issue of the variance, i.e. of how broad is the distribution of values for 5/5^°) around 
its mean value {S/S''-^'>). This dispersion is due to both the fluctuation on the averaging surface and to those due to 
ensemble fluctuations. Let us thus define: 



Var[5/5(°)] = ^(^/^W - (5/5(o)))'^ = {{S/S(o)y) - ((5/5(o)))' . (4.9) 
Inserting the definition (4.2) we get, after a little algebra: 

Var[5/5(°)] = (^^(1 + /.))((1 + m))-i - [{a{l + + {f,))-^) ' . (4.10) 

If we now make the same assumptions as before on expanding /i and cr, we find that the second term in (4.10) is at 
least of fourth order. Therefore, for the leading term in the variance we find the amazingly simple result (see also 
[20]): 



(4.11) 
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As in the case of IBR2 we only need to know the first-order perturbation, but this time the effect is completely 
independent of the averaging measure. As we shall see, the dispersion (which is actually the square root of the 
variance) turns out to be larger than the averaging corrections due to IBR2. 

We may note, at this point, that we could have also considered the dispersion of the angular average (5/S'(°-') due 
to the stochastic fluctuations, namely: 



Var'[5/5(°)] = ((5/5(0)) - (5/5("))] ^ ((5/5(o)))' - ((5/5(o))) , (4.12) 

which, after calculations similar to the ones above, gives 

Var' [5/5(") ] = (^. (4.13) 

Such a quantity is much smaller than the previous one, as can be inferred from Section VI, indicating that the main 
reason for the dispersion lies in the angular scatter of the data rather than in their stochastic distribution due to the 

ensemble. 

Let us finally identify the quantities appearing in IBR2 for our particular case. To this purpose let us stress that, to 
describe the real impact of inhomogeneties on the observational data, we should consider the scalar 5 corresponding 
to the true observed quantity. In the case of the supernovae data [.] this quantity should be the received flux of 
radiation, which is proportional to ~ d^^- On the other hand, if we consider the average of S = instead of the 
average of the flux, the difference is only of second order in our perturbative expansion ^. Since in this paper we will 
not evaluate such genuine second-order contribution to the average we can limit ourselves, for the sake of simplicity, 
to consider hereafter 5 = dL{zs,0°-). 

In such a case, from Eq. (3.31) we immediately find 

ai = Ai + A2 + A3 + Ai + A5, (4.14) 

where we have defined: 



A, = A2^2^^.: ^3=^l-^j/+; ^4 = -^l-^j/.; A, = -J2. (4.15) 

On the other hand, in our particular case, the general measure cP^ will be given by Eq. (2.9): 



d'^i = d'6^-i{wa, t(z„ ?«o, e-),e-) . (4.i6) 

We need to transform the integral in (P9 (over the "2-sphere" E) to the standard polar coordinates d'^9 = d6d<j) of 
the Newtonian gauge. To first order, it is easy to check (by taking into account the Jacobian determinant of the 
transformation 6'" — >■ 0", see Eq. (3.8)) that: 

d^9^ = J dcj)de sine i[asr,]{z„e'')f (1-2^,) . (4.17) 

The normalized unperturbed measure is then given by (dV)^°^ = d^n/An, where d^Q = d(j)d0 sin 9. Proceeding as in 
the previous section, and using in particular Eq. (3.29), we also obtain: 



Ml = -2*, + 4*av + 2 (^1 - Jizs.e-") = 2(Ai +A2+A^ + A4) . (4.18) 

The second-order induced backreaction IBR2, and the variance of d^/dj^^^ , can now be obtained by inserting 
the results (4.14), (4.18) into Eqs. (4.6) and (4.11), respectively. Their explicit evaluation, first analytic and then 
numerical, will be presented in the next two sections. 



^ Since the flux is proportional to dj^^{zs, 9°") the distance modulus (see Sec. VI) should be computed as a function of {{d^^)) ^1"^ rather 
than of (d^). Up to the second order the difference between the two is shown in the following expressions: 

= dl^^^ (l-l- -t- IBR^ - ^(^) ; (dZy = dl^^^ + + IBR,) . 
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V. INDUCED BACKREACTION AND DISPERSION: ANALYTIC CONSIDERATIONS 



This section will be devoted to a systematic evaluation of the various terms contributing to the induced backreaction 
IBR2 defined in the previous section, as well as to the variance associated to the (angle and ensemWe)-averaged value 
of 

Let us start by noting that the simplest way to implement the ensemble average of our stochastic background of 
scalar perturbations ^ is to consider their Fourier decomposition in the form: 



1 



(27r)3/2 



(5.1) 



where - assuming that the fluctuations are statistically homogeneous - _E is a unit random variable satisfying E* (k) 
E{—k) as well as the ensem&Ze-average condition: 



E(ki)E{k2)^5{ki + k2). 



(5.2) 



According to the results of Sec. IV, all corrections we need to compute are bilinear terms in the potential ^ 
always occurring in the combination {AiAj) or {Ai){Aj), where the quantities Ai, Aj are defined in Eq. (4.15). Let 
us illustrate a typical computation with one of the simplest contributions associated to Ai = —^s- Following the 
notations of Sect. IV we will denote with the angular brackets the integration over the two-surface S embedded on 
the light-cone (at w = wq, z = Zs). The measure of integration is the unperturbed normalized one, i.e. (Prt/Air. We 
then obtain: 



d^k d^k' 
(27r)3 



E{k)E{k') 



d^n 

in 



d^k 
(2^ 

d^k 
(2^ 
d^k d^k 



(27r)3 
d^k . 



E{k)E{k') 



2 

dk 
~k 

d^n 



(5.3) 



(27r)3 

d^k 
(27r)3 



An 

d(cos6>) .fc^^^Qgg 



1 



d'^il' 



sin(fcAr7) \ ' 
kAr] J 



°°dk /sin(fcA^y 



(5.4) 



where in the second line of both terms we made use of isotropy (i.e. 5'^ only dependent on A: = |fc|), and defined 
9 and 9' as the angles between k and x = rx and between k' and x' = rx' . We have also introduced the (so-called 
dimensionless) power spectrum of defined in general by: 



P^{k,r^)^^\^>,[ri)\\ 



(5.5) 



More complicated examples, that contain almost all the subtleties of these computations, will be presented in Appendix 
A. 

In general we have many contributions like the above ones, appearing in both the induced backreaction and the 
variance, and generated by all the Ai terms of Eq. (4.15). In this paper we will consider the particularly simple 
case of a CDM-dominated background geometry, with a time-independent spectral distribution of sub-horizon scalar 
perturbations, = 0. In such a case all the relevant contributions can be parameterized in the form 



{MA,) = 



dk 



dk 



P^{k)Cij{k,ria,T]s), 
P^{k)Ci{k,riQ,'i]s) Cj{k,rio,r]s 



(5.6) 
(5.7) 
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(where Cy , Ci are constant spectral coefficients) , valid for any given model of (time-independent) scalar perturbation 
spectrum. With such parametrization, the leading-order induced backreaction, called IBR2 in the expansion (4.6) of 
{(11) /df^^^'^ , can be written as 



4 5 



IBR2=2^^[(A,A,)-(A,;)(A,) 

The dispersion (from Eq.(4.11)) takes instead the form 



/ - P*(fc)2^^[c.,(fc, Vo,Vs) -Ci{k,i]o,T]s) Cj{k,riQ,r]s) 
■'0 1=1 ,=1 



Var 



^FLRW 



1/2 



= A/ (0-?) = 





1/2 


j=i j=i 





5 5 



Jo .-1 „_i 



1=1 j=l 



1/2 



(5.8) 



(5.9) 



We have analytically computed all the required coefficients Cij , Ci , and their final values are given in Table I for 
Cij(k,rif),rjs), and in Table II for Ci{k,r]Q,r]s). In Table I we also show the small-fc limit (kAri <^ 1) of C^j and of the 
products CiCj. For notational convenience we have introduced in the tables the dimensionless variable I — kAr], and 
we have used the definition: 

rl 



Sinlnt(/) 



dx 



■ smx. 



(5.10) 



TABLE I; The spectral coefficients Cij{k,riQ,rjs) for the (AiAj) terms defined by Eq. (5.6). We also give the kAr] <^ 1 limit 
(up to leading order in kA-q) of dj and of the products dCj for the coefficients defined in Table II. 



AiAj 



Cij{k,rio,Tjs) 



Cij for kA-q < 1 



Ci Cj for kAri < 1 



AiA, 



1 



AiA: 



-fSinInt(/) 



-2 



AiAi 



1 - 



sin(0 
I 



1 - 



■HsAi; ; 6 



- 1- 



-HsAr] I 6 



-ISL I ^] 11 

A?, I KsAr; J 3 



-f^ 1 ^11 

Ai; ^ «sAi; J 3 



AlAt 



Bin(i) 
I 



AoA: 



^ [-l + cosZ-f ZSinInt(0] 



4 - - 

9_ 



A2A3 







1 - 



KsAi; ; 3 



A2A. 



2 1 



fo+fs 

W3A?) J A77 



[1 



/0 + /-, 
Ar? 



H 

HsAri I 3 



A?) 



KsAj) 7 3 



3 



4+ (4 + ;^)cos/ + ZsinZ + /-^SinInt(Z) 







2 1 



1 

■HsAi7 



I 



WsAr, 



1 

HaAri 



A3A 



COS I 



sin(0 



Ar; 



1 



As A; 



sin(Z) 

i 



_ (1 ^ il 



A4A 



^^-^(2cos« + (-2 + Z 



Ari^ 



2 ^ sinj_ > 
I ) 



■H3A1; 



■H3A1; 



HsAri 



fo+3fs 
Ar) 



COS I + 



fo-f. Bin; , (/0 + /a)(-2 + iSinInt(i)) 



Ar) i 



+ 



Ar7 



fo-fs 
Ar? 



1 - 



HsAri I 3 



[-24 + 20r + (24 - 21^ + r) cos(/) 
+/(-6 + sin(0 + ?^SinInt(0] 



It should be noted that the computation of some of the above terms requires the explicit expression of the scale 
factor a(?7). In those cases we have assumed a dust-dominated scale factor with a(r/) = o,ivo){v/Vo)'^j ^-nd we have 
defined 



fo,s 



drj— 



V0,s Vin 



1 



(5.11) 



(recall that ?7i„ satisfies, by definition, r]in <C ?yo,s)- Using such a CDM-dominated model also imposes the following 
(zeroth-order) relation between Zs and Arj: 

2 



Ar] 



-1/2 



(5.12) 
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TABLE II: The spectral coefficients d for the (At) (Aj) terms defined by Eq. (5.7). 



A, 


Ci(fc,??0,T?s) 


Ai 


sin I 
I 


A2 


-fSinInt(0 


A:i 




Ai 




As 






This relation will be used in section VI when performing the numerical integration over fc as a function of Zg . 

In order to proceed with some qualitative and quantitative considerations it is important, at this point, to specify 
the properties of the power spectrum Pq,{k). Limiting ourselves to sub-horizon perturbations, and considering the 
standard CDM model, we can simply obtain ^'^ by applying an appropriate, time-independent transfer function to 
the primordial (inflationary) spectral distribution (see e.g. [ ]). The power spectrum of the Bardeen potential is 
then given by 



P*(fc) = ^^j A^T^ik) , Ai^A^-j , (5.13) 

where T{k) is a constant transfer function which takes into account the sub-horizon evolution of the modes re-entering 
during the radiation-dominated era, and is the primordial power spectrum of curvature perturbations, amplified 
by inflation, outside the horizon. The typical parameters of such a spectrum, namely the amplitude A, the spectral 
index and the scale fep, are determined by the results of recent WMAP observations ['^]. In our computations we 
will use, in particular, the following approximate values: 

^ = 2.45x10"^, =0.96, fco/ao = 0.002 Mpc"^ . (5.14) 

Finally, since we are mainly interested in the overall magnitude of the transfer function, it will be enough for our 
needs to approximate T(fc) with the effective shape of the transfer function for density perturbations^ with no baryon 
oscillations, namely T{k) — To(fc), where [2!)]: 

^°('^)= Lo+^°Co(.) ' ^"(^)=M2e + 1.8,), Co(.) ^ 14.2 + , .= ^, (5.15) 

and where fcoq is the scale corresponding to matter-radiation equality. We can easily see that the above transfer 
function goes to 1 for k ^ k^q, while it falls like k~^ log A: for ^ fccq- In the next section we will numerically evaluate 
the coefficients of Tables I and II, and we will discuss the corresponding backreaction effects on (cIl), together with 
its dispersion, for the scalar perturbation spectrum described by Eqs. (5.13) and (5.15). 



VI. INDUCED BACKREACTION AND DISPERSION: NUMERICAL RESULTS AND DISCUSSION 

In the following computations we will set ap — 1, — 1, and we will use h = i/o/(100kms^^Mpc^^) — 0.7. In 
that case we obtain [29] kgq ~ 0.036 Mpc~^, and we can more precisely define the asymptotic regimes of our transfer 
function as Tq ~ 1 for k < 10~^Mpc~^, and Tq ~ fc~^logfc for k > 2.5Mpc~^. This second scale is already deep 
inside the so-called non-linear regime, roughly corresponding to k >, 0.1 /iMpc^^ (see e.g. [30]) . 

Let us start by comparing the behaviour of the power spectrum with the behaviour of the coefficients in Table I and 
II. As shown in Table I, all the coefficients Cij and CiCj will go at most as 0(1) while their combination Cij —CiCj will 
vanish at least as 0{k^Ar]^) in the IR limit kAr] <^ 1. As a consequence, the infrared part of the terms in the induced 
backreaction and in the dispersion will give a subleading contribution, and we can safely fix our infrared cut-off to be 
k = Hq - i.e. we can limit ourselves to sub-horizon modes - without affecting the final result. 



The relation between density perturbations and metric perturbations is of course fully under control in the linear perturbative regime, 
but could be unreliable in the high-A: non-linear regime. 
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Furthermore, as can be seen in Table I and II, all the coefRcients Cij and CiCj will go at most as ©(fc^Ary'^) in the 
UV limit kArj ^ 1. As a consequence, all the integrals involved in the induced backreaction and in the dispersion 
will be UV-convergent and their main contribution will come from the range l/Arj <^ k < 2.5Mpc~^. In this range, 
some of the backreaction coefficients grow as positive power of k while the transfer function is not yet decreasing 
as log k/k^. In particular, there are only two leading contributions corresponding to the integrals controlled by the 
coefRcients C44 and C55 (this is the reason why the dispersion of the angular average of cIl, (4.13), is smaller than the 
one of d]^ itself (4.11)). For such coefficients we have, in fact, the following behaviour for kArj 3> 1: 



C44 
C55 



1 



k^AT]^ 



15 



SinInt(A:A?7) . 



(6.1) 
(6.2) 



From the above expressions it is easy to understand that (A4A4) will give the largest contribution for <C 1 (see 
Fig. 1), while will give the largest contribution for ^ 1 (see Fig. 2). 

It is also clear that both the induced backreaction IBR2 and the dispersion depend in principle on the UV cut-off kijv 
eventually used to evaluate the integrals. On the other hand, when kjjv is taken inside the regime where the spectrum 
goes like (logfc)^/fc^, the dependence on the cut-off will not be too strong. In particular the leading contribution to 
IBR2, given by 2(^4^4), depends very weakly on the particular value of kuv The leading contribution to the 
dispersion (5.9), controlled by (^15^5), has a somewhat stronger dependence on kuy (because of the extra power of k 
in the integrand of (A^A^)). The numerical integrations of {A4A4) and (^5^5) are presented in Figs. 1 and 2, where 
we illustrate the magnitude of the backreaction effects as a function of the redshift (in the range of values relevant 
to supernovae observations), and as a function of its dependence on the cut-off (ranging from kjjv — O.lMpc"^ to 
kuv — +00). We should emphasize that the two abovementioned contributions have a clear and distinct physical 
meaning in the Newtonian gauge. Going back to their explicit expressions it appears that (^4^4) represents a Doppler 
effect while (A5A5) is associated with lensing, both in qualitative agreement with previous claims in the literature 
[19, 21]. Note that the latter term can only appear in the variance, because it contributes to while lyA^A^i appears 




0.01 



0.02 



0.05 



0.10 



0.20 



0.50 



1.00 



2.00 



FIG. 1: The result of the numerical integration of Eq. (5.6) for {A^A^) is plotted as a function of for three different values 
of the UV cut-off: fe = 0.1 Mpc~^ (thick blue line), fc = 1 Mpc"'^ (thin purple line), k = -|-cx3 (dashed red line). Also shown (by 
the corresponding dotted curves) is the plot of the contribution to IBR2 (^44^4) — {Ai^^A^). 
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also in IBR2 sin ce A4 is pr esent in both /ii and ci. In Fig. 1 we have also plotted the full contri bution t o IBR2, 
namely (^4^44) — (^4) (^4) (see Eq. (5.8)): we can see that the difference from the behavior of the {A4A4) term is 
very small, with no qualitative effect on our discussion. 
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FIG. 2: The result of the numerical integration of Eq. (5.6) for {A5A5} is plotted as a function of Zs for three different values 
of the UV cut-off: k = 0.1 Mpc~^ (thick blue line), k = 1 Mpc~^ (thin purple line), k — +00 (dashed red line). 

It is important to stress that, although all considered backreaction contributions are (at any redshift) UV-finite, 
they can induce relatively big effects at the distance scales relevant for supernovae observations (see Figs. 1 and 2), 
provided the given spectrum is extrapolated to sufficiently large values of k. We also stress that the decoupling of 
the small-distance (i.e. high-fc) scales from the considered large-scale backreaction is due to the efficient suppression 
of the linear perturbation modes inside the horizon, an effect well described by the transfer function (5.15). 

Let us now sum up all contributions to the induced backreaction (5.8) and to the dispersion (5.9), and compare the 

results for (d^) ± d'£^'^ \J {af) with the homogeneous luminosity-distance of a pure CDM model and of a successful 

ACDM model. We will include into (di) only the IBR2 contribution (4.7), namely we wih set (di) = d£-°*^(l + IBR2). 
It is clear that a full computation should include additional contributions, of the same order as IBR2 or even larger, 
arising from second-order perturbations of (we are referring to the term called {(72) in Eq. (4.6)). Nonetheless, 
we believe that (modulo cancellations) our computation may estimate a reliable "lower limit" on the strength of the 
possible corrections to the luminosity-redshift relation in the context of our inhomogeneous geometry. 

The comparison between the homogeneous and inhomogeneous (averaged) values of d^ can be conveniently illus- 
trated by plotting the so-called distance modulus m — M = blogi^dL or, even better, by plotting the difference 
between the distance modulus of the considered model and that of a flat, linearly expanding Milne- type geometry, 
used as reference value (see e.g. [1!^]). In such a case we can plot, in particular, the following quantity: 



A(m - M) = 5 log 

The results are illustrated in Fig. 3 for the case of a cut-off kuv — O.lMpc"^, and in Fig. 4 for kuv — lMpc~^. 
The averaged luminosity-redshift relation of our inhomogeneous model is compared, in particular, with that of a pure 
CDM model and with that of a ACDM model with J7a = 0.1 and with JIa = 0.73. We have also explicitly shown the 

expected dispersion around the averaged result, by plotting the curves corresponding to {d]^)±d^^^'^ \J (af) (bounding 
the coloured areas appearing in the figures). 




{dL 



5 log 



10 



(2 + Zs)Zs 

2Ho 



(6.3) 
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FIG. 3: The distance-modulus difference of Eq. (6.3) is plotted for a pure CDM model (thin line), for a CDM model including 
the contribution of IBR2 (dashed blue line) plus/minus the dispersion (coloured region), and for a ACDM model with Qa = 0.73 
(thick line) and f^A = 0.1 (dashed-dot thick line). We have used for all backreaction integrals the cut-off k = 0.1 Mpc~^. 




FIG. 4: The distance-modulus difference of Eq. (6.3) is plotted for a pure CDM model (thin line), for a CDM model including 
the contribution of IBR2 (dashed blue line) plus/minus the dispersion (coloured region), and for a ACDM model with f^A = 0.73 
(thick line) and = 0.1 (dashed-dot thick line). We have used for all backreaction integrals the cut-off k = lMpc~^. 
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Let us stress again that the choice of the cut-ofF may affect (even if not dramatically) the final result when the 
values of kuv are varying in the range (0.1 - l)Mpc"\ while the precise choice becomes less important at higher 
values of kuy. We should recall, however, that the inhomogeneous model adopted in this paper is fully under 
control only in the linear perturbative regime, and that the spectrum cannot be extrapolated at scales higher than 
about k ^ IMpc^^ without taking into account the complicated effects of its non-linear dynamical evolution. The 
approximate coincidence between the above limiting value of kjjv and the scale marking the beginning of the non- 
linear perturbative regime is not only a particular consequence of the transfer function adopted in this paper (see Eq. 
(5.15)), but also an avoidable property of realistic non-linear perturbations spectra (see e.g. ["])■ 

As clearly shown by the two figures, the corrections induced by IBR2 on the l umino sity distance of a homogeneous 
CDM model, even taking into account the expected dispersion of values around (rfi), cannot be used to successfully 
simulate realistic dark-energy effects. In the figures we have also plotted, for illustrative purposes, an example of 
standard ACDM model with Oa = 0.1, which seems to be compatible (within the allowed region defined by the 
dispersion) with the prediction of our simple CDM-I-IBR2 model, at least for sufhciently high values of the cut-off 
scale (see in particular Fig. 4). However, we warn the reader that it would be wrong to conclude that our backreaction 
can mimic a fraction of dark energy of the order of f^A 0.1, because a similar conclusion might be reliably reached 
only by averaging inhomogeneities on a background which already includes a significant amount of dark energy from 
the beginning. 

Let us conclude this section with some important comments. First, as already stressed, a consistent second- 
order computation of the backreaction should include the contribution of (172). This requires a full treatment of 
metric perturbations and coordinate transformations up to second order (work is in progress on this point [ ]). 
Nevertheless, we can easily see from Eqs.(2.18) and (2.19) that (0-2) contains, among others, contributions of the type 
(^5^5) already computed in this paper. The behaviour of this term, in the asymptotic regime kAr] 3> 1, is very 
different from the behaviour of terms like {A4A4) which give the leading contribution to IBR2: the contribution of 
(^5^5), in particular, grows at large redshifts, as illustrated in Fig. 2. Using our results for (^5^5) - and barring 
cancellations - we can obtain, for instance, the following numerical estimates: with an UV cutoff kjjv — lMpc~^ we 
expect \{<J2)\ ^ 1.5 x 10^'^ at z = 1 and |((T2)| JJ, 4 x 10^'^ at z = 2. These expectations do not change much if we 
increase further the UV cutoff: they only go up by about a factor two even if we send kuv arbitrarily high (at least 
within our form of the power spectrum). 

These examples suggest that a full computation of (0-2) could strongly enhance the overall backreaction effects 
at large Zs, with respect to the effects due to IBR2 discussed previously. Hence a full second-order computation, 
possibly joined to a reliable estimate of contributions from the non-linear regime^, appears to be necessary before 
firm conclusions on the correct interpretation of the data can be drawn. Also, the different behaviour of the different 
backreaction contributions, at small Zs and large Zg, could represent an important signature to distinguish the effects 
due to averaged inhomogeneities from the more conventional dynamical effects of homogeneous dark energy sources. 

The second comment is that, although a reliable estimate of the full backreaction on the averaged luminosity 
distance requires a full second-order calculation, some suitable linear combinations of averages of different powers of 
dL only depend on the first-order quantity cti (defined by the expansion of du)- As an example, one can show that 
the following equality holds at second order for any value of the real parameter a: 

JXdiW^) - c^WW^ =l-a+ ^^^W)- (6.4) 

This quantity can be plotted for a given inhomogeneous model, and compared with its (deterministic) value in a 
ACDM model, for various values of a. The result is that the two models disagree for generic a, leading to the 
conclusion that realistic inhomogeneities added to CDM lead to a model that can be distinguished, in principle, from 
the conventional ACDM scenario. In practice, however, we only have a single quantity measured by the supernovae 
experiments (basically the received flux of radiation) , and one cannot exclude that the two models happen to give the 
same result for that particular observable. 

VII. CONCLUSIONS 

Let us briefly summarize the main ideas and results of this work. Using a previously introduced gauge invariant 
light-cone averaging prescription, as well as an adapted coordinate system, we were able to write down an exact 



The importance of such a non Unear regime was also recently underlined in [ ], following a different approach to describe the impact 
of inhomogeneities on the supernovae observations. 
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expression for the averaged luminosity distance ^^(z) (or any function of it) as a function of the redshift z. In 
principle such an expression can be used to study the effect of inhomogeneities even outside the domain of cosmological 
perturbation theory. In this paper we have only attempted a first study of the leading-order effect generated by 
stochastic perturbations of the type produced by inflation on top of a pure CDM model. 

We have then been able to separate two distinct contributions to the backreaction. The first, that we called 
induced backreaction (IBR), is sensitive to our specific averaging prescription and originates from correlations in the 
fluctuations of the luminosity distance and those in the integration measure. As a consequence, IBR can be computed 
to second order (giving what we have called IBR2) by using results from linear perturbation theory. The second 
contribution is insensitive to the averaging prescription but needs a full computation of second order perturbations in 
the luminosity distance itself, something that we are leaving to future work. Linear perturbation theory also allows 
to compute the expected variance on the (angular and ensemble) average of di^{z). 

Our final integrals for these quantities are nicely behaved both in the infrared (showing explicitly that perturbations 
on scales much larger than the source's distance do not contribute) and in the ultraviolet, provided short-scale 
perturbations behave roughly as in simple, realistic models for the matter power spectrum. This does not mean, 
however, that the final result is insensitive to the detailed structure of inhomogeneities in the non-linear (or even 
non-perturbative) regime. 

Although the power spectrum has an overall normalization of order 10^^ (from CMB data), it induces a backreaction 
which is strongly enhanced by phase-space factors of order (fc*/-ffo)'' (with p a power that can be as big as 2 or 3), 
where is some typical scale appearing in the power spectrum. These enhancement factors can bring the IBR effects 
many orders of magnitude higher. However, if we use a rather low momentum cutoff (by insisting on staying inside 
the linear regime), IBR2 effects turn out to be way too small to mimic an appreciable cosmological constant (see 
for instance Fig. 3). The total backreaction could instead be larger, but a reliable determination of its magnitude 
depends on being able to carry out the full second-order calculation mentioned above, as well as on having some 
knowledge about the non-perturbative (short-distance) contribution to d^. 

We have also noticed that the expected (angular) dispersion oi dL is quite large (see Figs. 3 and 4), particularly at 
small and at high Zg (with a minimum in the region of intermediate Zs, i.e. Zg = 0.3-0.6), something that one should 
be able to check quite precisely once more data become available. 

We note that the dominant contribution to the dispersion at small z comes from the {A4A4} term, and is therefore 
associated with "Dopplcr-type" contributions, since the quantity Ir is related to the source and observer velocities. 
This is in agreement with the claims made in [ ]. By contrast, the main contribution at large z comes from (Ac^A^), 
i.e. from terms usually referred to as "lensing" contributions. Again, this is in agreement with the literature (see 

The size and the nice IR and UV behaviour of our gauge invariant quantities can be contrasted with what happens 
if one approaches the backreaction problem using averaging on spacelike hypersurfaces, as in [25, 34]. In that case the 
calculation of the average expansion rate is UV convergent and practically cut-off independent. In a reasonable range 
of cut-off values it gives a backreaction of order 10~^ for the concordance model. However, for other quantities such 
as the variance of the expansion rate or the deceleration parameter, terms like ^ {d^'i')^ appear, giving UV-divergent 
integrals and a backreaction of 0(1) for relatively low cut-off. In our case, considering a UV cut-off equal to IMpc"^ 
we obtain a backreaction effect (estimated from (A5A5)) of the order of lO^'' for z ~ 1, namely 2 orders of magnitude 
bigger than the result presented in [25, -H]. Furthermore, there is no qualitative difference between the calculation 
of the induced backreaction and the one of the variance. They are all given by nice IR and UV-convergent integrals 
which never produce numbers of 0(1) even for very large momentum cutoffs. 

Finally, by suitably combining averages of different functions oi dL, we can again obtain results that, like the 
variance, only depend on linear perturbation theory. They show that, in principle, an inhomogeneous CDM model 
can be neatly distinguished from a homogeneous ACDM model. It remains to be seen whether any such combination 
is accessible to observations. Nonetheless, in the light of this last observation, it would appear quite unlikely that 
one will be able to fully account for the supernovae data in terms of inhomogeneity effects. Rather, depending on 
the full contribution of the second-order terms and, possibly, on the one from the highly non-linear regime, the effect 
of inhomogeneities could be relevant for future precise determinations of the critical fraction and equation of state of 
dark energy. In this spirit we plan to extend the calculations presented here to the (only slightly more complicated) 
case of generic ACDM models. 
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Appendix A. Other examples of backreaction integrals 



Here we present a detailed computation of some of the terms contributing to the induced backreaction and to the 
dispersion associated to the light-cone average of ^^(zs, ^''')- Such computations also define the corresponding spectral 
coefficients appearing in Table I and Table II. We will start with {A2A5) and {A2) (A5) as illustrative examples. We 
will then consider the two leading contributions given by {A4A4) and (A^A^) (and, for completeness, the associated 
terms {A4) {A4) and (A5) (A5)). 

Using equations (3.21), (3.26), (4.15), and working in the hypothesis of time-independent '^k, we obtain: 
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(A.2) 



In the second lines of Eqs. (A.l) and (A.2) we have applied Eq. (5.2) to remove the integration over fc', and we have 
used the isotropy of the scalar product k ■ x. 

We move now to the terms representing the two leading contributions to IBR2 and to the dispersion. Starting from 
Eq. (3.14), and using the same hypotheses as in the previous calculations, we have 
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Then, according to the definitions (4.15): 
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where 
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and where we have used the following exact integral result 
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In the same way we obtain: 
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Following a similar procedure Eq.(3.21) leads us to 
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Note that, for the same reason for which {A2) (A5) — 0, i.e. for the fact that the integrand of A5 gives zero when 
averaged over the angles, we finally also obtain (A^) (A5) = 0. 
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